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Abstract 

It is shown that the correlation functions of the random variables det(A — X), in 
which X is a real symmetric N x N random matrix, exhibit universal local statis- 
tics in the large N limit. The derivation relies on an exact dual representation of 
the problem: the fc-point functions are expressed in terms of finite integrals over 
(quaternionic) k x k matrices. However the control of the Dyson limit, in which the 
distance of the various parameters A's is of the order of the mean spacing, requires 
an integration over the symplectic group. It is shown that a generalization of the 
Itzykson-Zuber method holds for this problem, but contrary to the unitary case, the 
semi-classical result requires a finite number of corrections to be exact. We have also 
considered the problem of an external matrix source coupled to the random matrix, 
and obtain explicit integral formulae, which are useful for the analysis of the large N 
limit. 
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1 Introduction 



The spectrum of eigenvalues of complex Hamiltonians are often modelled by a random 
matrix theory, in which the random matrices belong to various ensembles according 
to the symmetries of the physical problem. The most common space-time symme- 
tries of the Hamiltonian lead to the consideration of ensembles of real, complex or 
quaternionic random matrices. In the simplest case one considers Gaussian prob- 
ability distributions. This simple choice is in many cases sufficient since it is now 
understood that the local statistics of the eigenvalues are universal, i.e. largely in- 
dependent of the probability distribution. The most commonly studied Gaussian 
ensembles , called GOE, GUE and GSE, are invariant under the orthogonal, unitary 
or symplectic groups, respectively, and they all have important applications fl], ||, 

In this article we follow our previous study for the GUE case of the characteristic 
polynomials of random matrices |Q. If X is an x random matrix , whose 
characteristic polynomial is det(A — X), we consider the average of products of such 
characteristic polynomials defined by 

k 

Ffc(Ai, A,) =< n det(Ai -X)>. (1) 
1=1 

In the GUE case we have derived in a previous article explicit formulae for those 
correlation functions, found then their asymptotic behavior for large N, and proved 
their universality in the short distance limit in which the differences Aj — Xj is the 
order of the mean spacing of the eigenvalues of X. As usual the orthogonal and 
symplectic cases are more difficult to handle. It turns out that there is a hidden 
duality in these problems between N the size of the matrices, and k the number of 
points in the correlation functions : we may turn the integrals over N x N matrices 
into integrals over k x k matrices and, since we are interested in large N -finite k 
limit, this is the required tool for obtaining the large N limit by the saddle-point 
method. We return below to the GUE case and exhibit its self-duality. However the 
GOE case turns out to be dual to the GSE. This duality, in the simpler case of all Aj's 

1 



equal, has been discussed recently within the orthogonal polynomial method 0, 
In both cases one may use geometry to reduce further the number of integrations. In 
the GUE case it relies on the Harish-Chandra-Itzykson-Zuber formula (HIZ) |[T^, |15 



For the GOE case it turns out that the N — k duality maps the problem into the 
GSE case, and the use of the geometry of the symmetric space U{2k)/Sp{k) leads 
to considerable simplifications. At the orders that we have considered one finds 
that there is a generalization of (HIZ) ; it is well-known that in that case the WKB 
approximation happens to be exact. For the GSE problem that one finds here, it is 
WKB, plus a finite number of corrections, which happens to be exact. Therefore we 
shall begin be re-exposing the unitary case at the light of this duality and of the HIZ 
formula. We shall then proceed to the GOE ensemble. 

2 Survey of the unitary ensemble 

For the Gaussian unitary ensemble (GUE), the random matrix X is a complex N x N 
Hermitian matrix, with a probability distribution function 

P(X) = lexp(-|trX2). (2) 

The average < ... > means integration with the normalized weight P{X) , with the 
Euclidean measure YiidXuYii^^j d^Xijd^Xij. It is easily shown that the F2{Xi,X2) 
reduces, up to a trivial factor, to the the kernel [§] Kj^{\i, A2) which characterizes the 
correlation functions of the eigenvalues of X. When all the X/s are nearby , i.e. in 
the short distance scaling region in which N is large and the products N{Xi — Xj) are 
finite, F2k{Xi, X2k) becomes, within an appropriate scaling, a universal function, 
i.e. independent of the specific distribution P{X). When all the Xj are equal, 

F2k{X) = F2k{X, A) =< [det(A - X)]^'^ > (3) 
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the 2k-th moment of the characteristic polynomiaL In the large N hmit, we have 
derived earher 0, |] 

F^k{\)=lk[2^Np{\)f, (4) 

in which p(A) is the density of eigenvalues, and 7^ is a universal factor. This number 
had been first computed for the circular unitary ensemble by Keating and Snaith 
who used the Selberg integral formula |p. 

There are several different derivations for those results. Let us here expose the 
duality which was mentioned in the introduction. We introduce Grassmann variables 
c and c, normalized to 

" dcdce'^'^" = I. (5) 
Then the characteristic polynomial may be written as 

N 

det(A - X) = / n rfc,cic„e'^^"(^^'^''-^'^^)^''. (6) 

•' a=l 

Repeating this k-times 

k . N k N k . 

The (normalized) integration over X , in presence of a matrix source Y, yields 

J rfXe-f = e-f . (8) 

We may now apply this to the matrix Yab = — J2t=i ^aaCba, generated by (^. Then 
one finds easily that 

k N N 

Tr(Y^) = - ^ Cb/jCfaa = -tr(7^) (9) 

af3=l a=l b=l 

with 

N 

lap = Y ^aaCa(3 ■ (10) 
a=l 

Our notation are as follows : "Tr" refers here to N-dimensional space, whereas "tr" 
refers to matrices acting in the k-dimensional space. We introduce next an auxiliary 



matrix k x k hermitian matrix B, such that 

gf tr72 = j dBexp (-ytrB^ + Ntr7B), (11) 

integrate over the Grassmann variables (which are now decoupled in the original 

N-dimensional space) and end up with 

r N 
Ffc(Ai---Afc) = J dB det{A - iB)^ exp -{—tiB^) (12) 

= e^''^' J dB{det B)^ exp -(y trB^ + iNtrAB) 

in which A is the diagonal matrix (Ai, ■ ■ ■ , A^). The problem is thus mapped into 
Gaussian integrals over k x k hermitian matrices as announced. This dual represen- 
tation is of course well adapted to the fc-fixed, A^-large, limit that we are considering, 
since (|1^) contains k"^ variables instead of A^^ in our starting point. It is not difficult 
to proceed from (12) and derive the scaling results which were given in 0]. 



However it turns out that it is simpler , and necessary in view of what comes next 
for the GOE problem, to integrate out the unitary degrees of freedom in (|12|). This 
is done through the HIZ formula |T^, IT3 which gives the integral over the unitary 



group U{k) : 

iTvtrUXUtY ^ ueLi<jj<fce 

"f^e = Oat— — — [id) 

A{xi, - ■ ■ ,Xk)A{yi, - ■ ■ ,yk) 

in which the Xj's and y/s are the eigenvalues of the Hermitian X and Y respectively; 
A(a;i, ■ ■ ■ , Xk) is the Van der Monde determinant 

A{xi,- ■ ■ ,Xk) = Y[{xi - Xj). (14) 

i<j 

It is well-known that the formula (0) happens to be exact semi-classically, i.e. if 
one retains only the sum of the k\ stationary points in the space of unitary matrices, 
weighted by the Gaussian fluctuations around each of them. Higher corrections hap- 
pen to cancel exactly. This leads immediately to an integral over the k eigenvalues 
bj of B, rather than over the fc^ matrix elements : 



F,(A„...,A.)=cet-^^ n^^^^r e--w^^-^^^'^on^^^ (is) 

\Z=1 / Kl' - ^l') 



When we consider simply the k-th moment of the characteristic polynomials, namely 
the case in which (A = Ai = -- - = Afc), the previous formula reduces to 



F,(A) = F,{X, . . . , A) = Ce^'-' j [X{dh,h^ ^^.(^f-^^^^A) ^^^^ _ 

\l=l ) Kl' 



(16) 



Note that the above representations of the correlation functions of the characteristic 
polynomials, in terms of integrals over k variables, are exact for any size N x N of 
the random matrices. It is then simple to find the large-N limit of those functions by 
saddle-point integration. If we focus to even values of k, and substitute 2k to k (the 
odd case is doable of course, but it leads to an oscillatory behavior) the saddle point 
equation for each 6/ is 

- iXih -1 = (17) 

whose roots are hf = (iXi ± — \f)/2. In the scaling limit, in which the A; — A;' 
are of order 1/A^, the eigenvalues bj are very close and one must pay attention to the 
Vandermonde determinant in the integration measure. Finally the leading saddle- 
points correspond to equal numbers of bi close to either b~^ or b~ , with 6^ = {i\ ± 
-\/4 — A^)/2. There are thus f = 2k\/k\k\ saddle-points of equal weight. The 



. k , 

combinatorial factor 7^. of (H) is then simply 



7fc 



f2k\ {h 



Ik 



\2 



\k J h 

k-l 

n 



'■2k 

l\ 



1=0 



{k + l)\ 



k 
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where we have used 

k 

= (19) 
1=0 

(This formula is used when we consider the gaussian fluctuations near the saddle- 
point in which k of the 6j's are near b^ and the other half are close to the saddle 
point b~ .) 



Finally this representation through an integral over a finite matrix matrix B, 
may be generalized to the case of an external matrix source A coupled to the random 



matrix X B 11, 12 



2 trX^+NtrAX 



f dxf[det{Xi-X)e- 
1=1 

1 N 2k 2k 

= l^JdkUm^^-bi)mbi-bi')e-^^'^^'''^''^ (20) 
We shall now transpose these techniques to real symmetric random matrices. 

3 Real symmetric matrices and characteristic poly- 
nomials 



Again we consider 



k 

Fk{Xu Afc) = J rfXe-T*"^^' n det(Ai - X), (21) 



i=l 



in which X is a real symmetric N x N matrix, X = X . It is worth remembering that 
real symmetric matrices form the Lie algebra of the symmetric space U{N)/0{N) 
(the Lie algebra of U{N) consists of A^^ complex hermitian matrices ; the imaginary 
part of those matrices are the A^(A^ — l)/2 antisymmetric generators of 0{N) ; the 
real parts are the A^(A^ + l)/2 real symmetric generators of the coset) 

Using again Grassmann variables, and the representation (|^) of the characteristic 
determinants we are led again to an integration over real symmetric matrices in the 
presence of the matrix source Y = — Y.t=i CaaCba- This gives 

j ^Xe-^TrX^+iNtrXY ^ g-f Tr(Y2+YYT) _ ^22) 

We have dealt earlier with 

Tr(Y2) = -tr(72) (23) 

which led to the integral ([TT| ) over an hermitian k x k matrix. In addition we have 
here 

Tr(YY^) = tr(UV) , (24) 
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with the matrices U and V defined by 



N 

Ual3 = CaaCap (25) 
a=l 



N 



K/3 = X] CaaCa/3 (26) 



a=l 



Defining the complex conjugation of Grassmann variables as (C1C2)* = C2C1, we 
have 7 = 7^, = [/. Therefore, we may again decompose the remaining quartic 
terms in the c's and c's as 

g-f tr(UV) ^ y ^2^g-mr(DtD+ivtD+iDtV) (^^7) 

where D is a complex k x k antisymmetrix matrix, D = — D^. Then we have 

X J (^i^(^2:)e-''Vtr(B2+DtD)+N^CaaCa/3B/3a 

Those auxiliary matrices B and D allow us to integrate over each pair Ca, Ca indepen- 
dently of the other pairs. It is convenient to define 

^a=f?V (29) 



Ca. 

For a given antisymmetric matrix M, (M = —M'^), we then have the following 
formula : 

iN 



= [-PfM]^ = [detM]^/2 (30) 

where Pf is the pfaffian of the antisymmetrix matrix M. Applying this to our prob- 
lem, we deal here with 

*^=(-(A-iB) Dt ) (31) 



in which A = diag(Ai, ...Ak). Thus we finally obtain 

Ffc(Ai,...,Afc) = J rf5c/De-^*'-(^'+°'°H-Pf^]'^ (32) 

This representation ( p2D of Fk in terms of a finite number of integrals, here 2k'^ — k 
integrals (one k x k hermitian matrix B, one complex antisymmetric k x k matrix 
D), again exact for any A^, is a solution to the problem. However, contrary to the 
GUE case, it turns out out that a direct use of the saddle-point equations fail in the 
scaling limit. In other words every term of the perturbative expansion around the 
saddle-point turns out to be relevant in the regime in which the products iV(Aj — Xj) 
are finite. 

One may use the remaining invariances of this representation, to reduce further 
the number of integrations. A unitary transformation U, among the c^, which diag- 
onalizes the Hermitian matrix B, one of the block matrices of M, transforms D into 
D U*DU'^] in other words one can diagonalize B and keep for D an antisymmetric 



matrix. Therefore, applying the Harish-Chandra-Itzykson-Zuber formula [15| for 
the integration over the unitary group, i.e. over the relative unitary transformation 
between the diagonal matrix A and the eigenbasis of B, we obtain 



1=1 



X (-PfM)^ (33) 

where now the matrix i? in M is diagonal ; we have reduced the integrations to fc^ 
variables, instead of 2k^ — k. However it turns out that this is still unsufficient : the 
inapplicability of the saddle-point method in the scaling limit is still a problem if we 
proceed from (^). It is thus necessary to return to the underlying geometry of the 
space of matrices M in the representation (0). In order to make the quaternionic 



structure more apparent we return to (|28D and define the spinor 



^aa =[;"] (34) 
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and the adjoint 



(35) 



Then the quadratic form in the Grassmann variables of (pS]) takes the form (repeated 
N times for each index a that we drop) 

a, (3=1 a, (3=1 

(36) 

in which the qa(3 are quaternionic matrix elements, i.e. linear combination of the 
Pauli matrices. The identification in terms of 2 x 2 matrices is thus 

Ba,(3 ~iD*, 



qa,(3 



-iD^,p El 



a,(3 
a,(3 



This defines a self-dual quaternion matrix |jT], || , i.e. 

qi(3 = 5'/3" 



(37) 



(3^ 



and qa(3q(3a is a multiple of identity. Let Mq be the quaternionic matrix whose elements 
are the quaternion qa(3 and M the quaternionic matrix with elements. 



Ma(3 — qa(3 + i^^a^aP- 



(39) 



The Grassmannian integration leads to 



Yl dcadca exp -— J2 ^» (^a,(3 + iK5a(3) ^(3 = Q det M = -PfM, (40) 



a,(3=l 



in which Q det denotes the quaternionic determinant |[T|, §| . In addition 

tr{B^ + DD^) = tr(M^) = 5]q„^q^„. 

a(3 



(41) 



It may be clarifying to show this quaternionic construction explicitely for the 
k = 2 case. There one has 



M 



( ° 


d 


Ai - iBn 


-iB2i 


-d 





-iBi2 


A2 — iB 


-Xi+iBu 







-d* 


V iB2i 


-X2 + iB22 


d* 






22 



(42) 
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and 

- PfM = |dp + (Ai - iBn)(A2 - iBas) + |B 

The equivalent quaternionic construction is 

^ ^ / gii + iXi qi2 \ 
\ q2i q22 + i>^2 ) 

with 

gn = Snl , qi2 = {^B^2)l + ^{'^B^2y3 + i{^d)a^ + t{'^d)(X2 

q2i = qi2 , 922 = ^22! (45) 

and 

Qdet{M) = (gii + i\i){q22 + i>^2) - ququ 
= {Bn + ^Xl){B22 + 1X2) + \Bi2\^ + \d\^ = -PfM. (46) 

Therefore we end up for the correlation functions with the following duality : 

Ffc(Ai,...,Afe) = J rfM(QdetM)^e-^*'^(^«'). (47) 

The original integral over real symmetric N x N matrices is replaced by integrals 
over quaternionic matrices which depend upon (2/c^ — k) degrees of freedom. Those 
matrices are the generators for the symmetric space U{2k)/ Sp{k). 

This representation (|47|) in terms of a finite number of integration variables is 
a priori well adapted to the large N-limit. However it turns out that in the scaling 
limit of interest, the contributions of the non-gaussian fluctuations around the saddle- 
points are all relevant. Therefore it is necessary to eliminate first the "angular" 
degrees of freedom. When all the Aj are equal, namely if we consider the moments 
of the characteristic polynomials, one can simply diagonalize the symplectic matrices 
in terms of k eigenvalues and then proceed to the large N limit. This is done in the 
next section. However if the Aj's are unequal we need some equivalent of the HIZ 
formalism, which will be described afterwards. 
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(43) 
(44) 



4 Moments of the characteristic polynomials 

We first note the trivial k = 1 case : Fi{X) = (det(A — X)) is simply 



which, up to a trivial factor, is the Hermite polynomial Hn(V^A) which has an 
oscillatory behavior for large N when A belongs to the support of Wigner's semi- 
circle. Therefore we consider from now on the more interesting even correlation 
functions. 

When all the Aj's are equal, the matrix A is a multiple of identity, and one can 
diagonalize the quaternionic matrix M through a transformation belonging to the 
symplectic group Sp{2k). The transformation of M into the diagonal matrix T = 
diag(ti, ...tk) yields the Jacobian J = [A(t)]^, (A(t) is the Vandermonde determinant 




(48) 



This gives simply 



F,,{\) = (det(A-X)2'=) 

2k 2k 

= C I \{dU{X{urX{{U-tvYe-''T.^'^^^'^^T.^^ (49) 




1=1 1=1 KV 



The integral representation (|49D is well suited to the study of the large N limit. 

2k 

Exponentiating tf term as e^'°s*', the integrand is of the form exp(— X^ fi^i)) with 



1 




The saddle points for every ti are solutions of /'(tc) = 0,i.e. 



2t^ - 2iXt -1=0 



(51) 



The two solutions are given by 



^{tX±V2^) 



(52) 
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The difference (t^ — t ) is proportional to the semi-circular density of eigenvalues of 
the GOE ensemble : 

t+ - r = 7rp(A) (53) 



where p(A) = \/'2 — A^/vr. Expanding ti around either t+ or t^, we find that the 
leading saddle-points are those in which half of the ti{l = l,...,2k) are near 
and the remaining half near . (Other choices give oscillatory contributions in 

2k 

exp{—N fjti)) , which damp the large A^-limit). Therefore we have to add the 

1 

{2k)\/{k\k\) leading saddle-points corresponding to the distribution of half of the t;'s 
near t~^, and the other half near t~ . The measure term given by the 4-th power of 
the Vandermonde determinant, yields a factor (7rp(A))^'^^ from the k variables near 
and the k near . The exponent / is then expanded around t+ or , and the 
remaining integral factorizes into an integration around and another one around 
t~ . The integration around if:"'' is 

n dUe~^r'it-Kt-t-? fl^t^ _ t^y = ( ] r-'^ n(2/)!. (54) 

'1=1 ^<j jNf"{t+) 1=1 



Noting t+t- = -1, and t+ - r = np{\), we find f"{t+)f"{t-) = inp{\)f. We need 
to fix the normalization constant C in (l49|) . It is obtained from the integral, 

2k 1 2k 

I n^^^e-f i:*? X{{U-t^f = (^)^^'-'^n(20! (55) 

i<j 1=1 

The constant C in ( P?D is thus the inverse of this number. This constant appears 
as a normalization for the n-point correlation function of the Gaussian symplectic 
ensemble 0. Thus including this normalization constant, F2k{)^) becomes in the 
large N limit as 

Fi^''''^=^,N''\2npiX)r'^' (56) 



(2fc)![nti(20P^ 

k\k\ n'=i(20' 

A (2/-1)! 
fj^ {2k + 2l-l)\ 



7fc 



A (2^-1)' (.7, 
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For example, in the case of F2(A) =< (det(A — X)y >, it gives 

F2(A)(^^^) = l:N\2np{X))' (58) 

This result agrees with the result which one would deduce from lim p(Ai, A2) / (Ai — A2), 
where the limit means Ai A and A2 — > A, ( details are given in appendix A). 

The value of 7^ agrees with the result of COE (circular orthogonal ensemble) 
found by Keating and Snaith [^j through the Selberg integral formula. In the COE 
case, however, the density of state is a constant, and the factor p(A) is absent. This 
result is to be compared with the earlier result for the GUE, 

F^^^ = lk{2N7rp{X)f (59) 

where 7^ is given by ([T^) . This universal constant 7^ appears also in the average of 
the moments of the Riemann ^-function, and it has a number theoretical meaning 

5 Correlations of characteristic polynomials 



The integral representation ([T^ ) of the correlations functions -Ffc(Ai, A^) is not 
unitary invariant unless the Aj's are all equal. Therefore if we parametrize the matrix 
B as B = WbU, in which 6 is a diagonal matrix, we have to consider the HIZ integral 

/ dU exp ^trUtbUA = "^""^ (60) 
J A(6)A(A) ^ ^ 

which is well-known to be WKB exact . This explains why , in the GUE case, it is 
equally possible, to apply the saddle-point method with or without integrating out 
the unitary group. 

In the symplectic case we are not aware of any similar explicit result ; however 
it will be shown now, at least for the lowest values of k, that the integral over the 
symplectic group can be performed exactly. The result is remarkably that, in this 
case, WKB plus a finite number of corrections is exact. 
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In the k = 2 case, we have 

F2(Ai, A2) =< det(Ai - X) det(A2 - X) > (61) 
As shown in the previous section, it is given by the integral (551). We first evaluate 



explicitly the angular integral. We first diagonalize i? by a unitary transformation, 
and then write the eigenvalues bi and 62 in terms of new parameters t and c, 

h = (l-c)ti + ct2 

62 = cti + (l-c)t2 (62) 

Then we have 61 + 62 = + h, and 61 — 62 = (1 — 2c) (ti — t2). Since the integrand 
is a function of we change variable to = c(l — c)(ti — t2)'^ ; then we have 
bl + bl + 2|(ip = + tl, and 6162 — MP = tih- Finally since the parameter c is 
restricted to the interval < c < 1, we replace it by c = sin^ 9. 
This leads to 

r-l roc (+_ _ +„^3 



i^2(Ai,A2) = /' dc r dhdt2{l - 2c) W iht^f 
Jo J-00 [Ai — A2) 

X ^-N(tl+q)+2iN{tiXi+t2X2)-2iNc{ti-t2){Xi-X2) (^g3^ 



The integration over c yields 

F2(Ai,A2) = r dtidt2(tlt2)'^e-^(*?+*')+2^7V(t,A,+t,A,) 
J —00 

1 tl — ^2 n2 , ^ tl — h 1 ,N 

When Ai = A2 = A, it reduces as expected to (|^). This formula may be easily 
checked for finite values of , since it reduces to Gaussian integrals over ti and 
t2 ; for instance in the simplest case iV = 1, it gives -F2(Ai,A2) ~ A1A2 + 1, which 
agrees with the direct calculation (in this case the trivial integral over the real axis 
J dx{Xi - x)(A2 - x)e~5^^). 

However the representation (|6^), which is exact for any A^, makes it clear 
i) that the large A^-limit may be found through a saddle-point integration over ti and 
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t2 ; this will be done below. 

ii) that in the universal local limit of interest, in which goes to infinity, Ai — A2 goes 
to zero and A^(Ai — A2) remains finite, the large iV-limit could not have been taken 
earlier. If, for instance, we had used the saddle-point method at the level of (p3|), we 
would have missed the second term in the bracket of (|6^) . If, at the early level of (33), 



we had recognized that the regime of interest requires to expand beyond the Gaussian 
approximation to the saddle-point, it would have appeared unexpectedly that the 
expansion stops after the first correction. Therefore (|6^ ) could have been obtained 
by a semi-classical approximation with a finite number of corrections, here just one. 
This is analogous, although not as simple, to the Harish-Chandra-Itzykson-Zuber 
formula for the GUE case, which is semi-classically exact, without any correction 



term |T^. The above integration over c is therefore , for = 2, the corresponding 
HIZ formula for the symplectic group. 

For higher values of k we need a more elaborate strategy. The HIZ formula may 



be easily derived by considering the Laplacian operator ||16| , 



L = -IHFT- (65) 



Its eigenfunctions are plane waves 

One can construct a unitary invariant eigenfunction of L, for the same energy A^^trA^, 
by the superposition 

/ = |c^[/e^^™u^ (67) 

which is nothing but the HIZ integral. The integral beeing unitary invariant, it is 
a function of the k eigenvalues of X. The same considerations hold for the three 
ensembles (3 = 1,2 and 4, corresponding to the orthogonal, unitary and symplectic 
ensemble, with 

/ = / e''''''^''^~'dg. (68) 
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The Laplacian, expressed in terms of a differential operator on the eigenvalues ti 
reads 

with the eigenvalue e 

e^AT^EA^ (70) 
The t-dependent eigenfunctions of this Schrodinger operator have a scalar product 
given by the measure 

{'pM^ I dh---dtk\A(h---tk)f ipl{h---tk)^2{ti---tk) (71) 

The measure becomes trivial if one multiplies the wave function by |A|^/^ . Thus if 
one changes I{t) to 

iP{h---h) = \A{h---tk)f/'l{h---tk), (72) 
one obtains the Hamiltonian, 

E S - - 1) E 77^]^ - (73) 

i=i ^^i ^ i<:j ^j) 

For P — 2, the solution is again given by plane waves in the ti and (taking into 
account the symmetry under permutations of /), one obtains the HIZ formula. 

In the P — A case, the problem is less trivial, but simple for finite values of k. For 
k = 2, a. solution of this equation is 

^° + ;v(«.-«')(A.- A.) ' 

The symmetry of / under permutation of thed tj's leads then to the solution 



N{ti - t2) (Ai -X^y ' N{h - h) (Ai - A2) ' 

(75) 

Then, after multiplication by the Vandermonde factor, we obtain the required sym- 
plectic HIZ formula (for k — 2), 

1 
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ij (76) 



For general k {f3 = 4), the solution of (|73|) is of the form 
where x satifies 

E| + 2»ivi:A.|-E^lx = (78) 

1 = 1 I 1=1 ' 1<J V * J/ 

k Q2 4 

The operator ^ — ^ — ^ — ■ annihilates the function A^^(ti, ■ ■ ■ ,tk) . Con- 

i=l ^^i i<j V-i ~ tj) 

sequently the solution of (^) may be written 

in which f{ti ■ ■ - tk) is a polynomial of degree k{k — l)/2 in the t^s. Defining 

Tij = N{X,-Xj){ti-tj) (80) 

one finds for k = 3 



2 1 1 1 

x=[i --(— + — + —) 

i ri2 Tsi 



'T12T23 T23r3i r3iri2 ri2r23r3i 
Again, as for k = 2, one sees that the successive terms in the r.h.s. of (|81|) are of 
same order in the limit of interest, and again they could have been obtained through 
a finite number of corrections to a semi-classical calculation. It is remarkable that the 
series of x stops at the order the inverse of the Vandermonde ; thus the symplectic 
HIZ integral is expressed as the sum of a finite number of terms. The successive 
coefficients of each term are determined by the equation ([731). 

Using this modified HIZ formula for the symplectic case, we obtain for the k=3 
case, -F3(Ai, A2, A3) which is expressed by 

F3(Al,A2,A3) = / rftirft2rft3e"^(*?+*^+*3)+2.A^(Ai*i+A...+A3*3)(i^^^i3)^(|W )2 
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1 1 1 

+ 'U(Ai - A2)(ti - h) ^ N{\2 - A3)(t2 - ta) ^ iV(A3 - Ai)(t3 - hY 
1 1 



N^iXi - A2)(A2 - A3)(ti - h)it2 - ts) iV2(A2 - A3)(A3 - Ai)(t2 - - ti) 
' -4 .r..L.u. ] (82) 



iV2(A3 - Ai)(Ai - A2)(t3 - ti){h - t2) 2 iV3A(A)A(t) 
where A(t) = (ti — t2)(^2 — ^3)(^3 — ^i)- (Note that the symmetries of the integrand 
allowed us to keep only the single solution (^TJ), without adding permutations). For 
low values of (i.e. N=l or 2, one verifies easily this result by a direct integration 
over 1 X 1 or 2 X 2 matrices. 

Higher values of k may be handled in a similar way, but the combinatorics become 
quite heavy. For instance in an appendix the solution of the case A; = 4 is given 
explicitly and, although again it consists of a finite number of terms, it is quite 
cumbersome. 

As is now clear , those integral representations make it easy to find the scaling 
limit (large N, finite A^(Aj — Xj))- For instance for k = 2 one finds the saddle point 
values of ti and t2 from (p^) in the large N limit. 



U='-\. + y2-\^ (83) 

where i = 1,2 and e = ±1. We use the parametrization 

Aj = \/2 cos 9i. There are a priori four saddle-points given by (|83D, but the two 
dominant ones are ti = t2 = ^j^*^^*^ ^ ~ short distance limit, 

large and A^(^i — 62) = N612 finite , we obtain 

where 9 = {9i + 62) /2. The semi-circle density of states is given by p(A) = ^sin^, 
and 6*12 = — ts^^^-^- Thus we obtain in the scaling short distant limit, 

F,(VA,) = C|^-2i!Hl (86) 

where x = itN{Xi — A2)p(A). It is interesting to note that this function may be 

expressed as a half-integer Bessel function, since ( — —) = — J — r 73/2(0;). 

X'^ V 2x'^ 
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In the unitary case, the sine kernel is similarly a half integer Bessel function since 

sin X nr 



6 Extension to an external matrix source 

In the GUE case, when an external matrix source A is coupled to an Hermitian 
random matrix X, as we have discussed earlier, F2fc(Ai, A2fc) is given by (^0|) . The 
degrees of freedom provided by the eigenvalues of A are useful to study a number 



external source matrix A, we can study the problem of a closing gap in the spectrum 
of random hermitian matrices p. Thus it is interesting to consider this external 
source problem for real symmetric matrices as well . 

One can always assume that the external source matrix A is diagonal. In the 
method of integration over Grassmann variables used in section three, it is simple to 
include the external matrix A : 



where Y = — Y,a CaaCba- Since A is diagonal, the term trAY gives simply the extra 
term exp[iN J2 O'nCnaCna] in the integrand of in (pHf), . Therefore, repeating the 
calculations of section 2, we find that eg. (pOD is modified as follows : 



of new universality classes For instance by tuning the eigenvalues of the 




= e 2 



f tr[{Y-iA)2+(Y-iA)(YT-iA) 



(86) 




n[-PfM(^')] = n[detM(^')]/2 



(87) 



j j 



where PfM*-"-* is the pfaffian of the antisymmetrix matrix M*^"^ given by 




8. 



.8 
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in which A = diag(Ai, ...Ak) and ajl = diag{aj, aj). Thus we finally obtain 



N 



Fk{Xi,....,Xk) = J rf^rfZ^e-^^'^^^'+D'D) n[-PfM(^)] (89) 

i=i 

This integral can be expressed in terms of a quarternion matrix Q, which can be 
diagonalized by the symplectic group Sp{k). When all the Aj's are equal to a single 

A, we get 

k N k 

F,{\, A) = / n n(^' - n(^' - 1.) v^i:*?^^^^^!:*' n du . (90) 

1=1 j=l Kl' 1=1 

For the case of different Aj's , this formula is modified by an extra factor as in the 
previous section. 

As an example of the usefulness of the above representation, we choose an external 
source with only two opposite eigenvalues ±c, with half of the eigenvalues equal to 
+c and the other half to — c. This gives a factor (tf + c^)^/^ = exp ^ log (t^ + c^) in 
the integrand. Expanding it in powers of t^, the total coefficient of t"^ in the exponent 
vanishes for c = 1/^/2. Therefore in the large N limit, we obtain at this new critical 
point 

< [det(X)]'= >= J e-'^'^'^'dQ (91) 

where Q is a kxk symmetric quaternionic matrix, and the the average < ... > is eval- 
uated with the distribution in the presence of the external source whose eigenvalues 
are ±c = ±1 / v^. 

One could make other choices for the eigenvalues of the external source matrix 
A, and obtain thereby higher multicritical points with terms such as trQ'^ in the 
exponent, in analogy with the GUE case in an external matrix source |0 . 



7 Summary 

In this article, an exact representation of the k-point functions (na=i det(Aa — X)), 
averaged over N x N real symmetric random matrices, has been derived in terms of 
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an integral over quaternionic k x k matrices, invariant under the unitary symplec- 
tic group. This representation leads to an easy calculation of the moments of the 
characteristic polynomials (Ai = ■ ■ ■ = A^). In the large A^-limit one finds 

^-°'''^npi^^A'-'(2MA))--. (92) 

to be compared to the earlier result for the GUE, 

^ir''^ = n(^(2iV-P(A))'^^ (93) 

For unequal A^'s, in spite of the fact that the integral representation involves a fi- 
nite number of variables, in the large A^-limit the corrections to the saddle-point, 
in the scaling regime A^(Aj — Aj) finite, are not negligible. A generalization of the 
HarishChandra-Itzyson-Zuber formula is shown to solve the problem. Remarkably 
this formula is" nearly" semi-classical, in the sense that it happens that the semi- 
classical expansion terminates after a few terms, a number of terms which increases 
with k but not with A^. Then the saddle-point method may easily be applied for 
large N, and this leads to explicit asymptotic formulae for the correlation functions 
of the characteristic polynomials. Finally this may be generalized to include an ex- 
ternal matrix source in the probability measure. Real symmetric random matrices 
appear as models of numerous physical time-reversal invariant Hamiltonians. For 
instance the orthogonal matrix model with an external source has been investigated 



as a model of glassy behavior [0. The results of the present work for the moments 
and for the correlation functions in an external source may be of interest for such 
problems. 
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Appendix A: The solution for k = A 



From ( |78| ) and ([79D, the polynomial / satisfies 

k f)2 ^ Bf (91 ^ Bf 

,A.l, 

The solution of this equation is obtained by a perturbation expansion in powers of 
the Aj's, but it ends at the level of the Vandermonde A(Ai ■ ■ ■ A4). Using the notation 
of (pnp, Tij = N{ti — tj){\i — we obtain 

f = C[l- -{Ti2 + Ti3 + Tu + T23 + T24 + T3i) 

- ........ ...... 

+ + T14T24 + T24T34 + T23T34 + T13T34 + T13T23) 

- — (T12T34 + T13T24 + T14T23) 
io 

+ ■^{ti2Ti3Tu + n2T23T24 + T13T23T34 + TUT24T34:) 

6b 

+ T14T34T23 + TUT 24r 23 + r 12T24T 34 + ^ 12^-23^ 34 
+ T12T14T23 + Tl3''"l4T23 + Tl2Tl3''"34 + T12T14T34 
+ T-13^34T24 + T13T 24T 23 + T 14^ 24^13 + ^13^ 12^ 24) 

+ -^^iT 12T 23T"34T- 14 + T uT 13T 2aT 34 + ^13^" 14T" 24T 23 + ^-12^ 14T247' 34 + ^ 12^ IAT 24^ 23 
+ T12T I4T 24T 13 + Tl2T" 13^" 23T34 + Tl27" 13T-23T 14 + T 12^ I3T23T 24 + T I2T24T 23^ 34 
+ Tl4T24T23''"34 + T 13^24^ 23^34 + ^12^ I4T 13^ 34 + T 14^ 13^23^ 34 + T I4T I3T 34T 24) 

- 13T- 24 T" 23^" 34 + ^ 12^ uT 24^ I3r 23 + 1" I2r uT 24!" I3r34 
+ Tl4:Tl3T24,T23T34 + T12T14T13T34T23 + Ti 2X14X24X23X34) 

- 7^^^-12X13X14X23X24^34] ■ (A.2) 
Zoo 

The HIZ integral is obtained by requiring the symmetry under permutation of the 
tj's in the final expression for J, thus from ([75| ) (and a replacement of A by 2A), 

I = Cj^^^^^^ifih,...,t,)+perm.of f) (A.3) 
22 



where the last term means that one adds the terms in which one permutes the tiS 
for fixed A's . 
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